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ABSTRACT 

For the first time the cross-correlation between type la supernova host galaxies and surrounding field 
galaxies is measured using the Supernova Legacy Survey sample. Over the z—0.2 to 0.9 redshift range 
we find that supernova hosts are correlated an average of 60% more strongly than similarly selected 
field galaxies over the 3— 100" range and about a factor of 3 more strongly below 10". The correlation 
errors are empirically established with a jackknife analysis of the four SNLS fields. The hosts are more 
correlated than the field at a significance of 99% in the fitted amplitude and slope, with the point-by- 
point difference of the two correlation functions having a reduced for 8 degrees of freedom of 4.3, 
which has a probability of random occurrence of less than 3 x 10~^. The correlation angle is 1.5±0.5", 
which deprojects to a fixed co-moving correlation length of approximately 6.5 ± 2h~^ Mpc. Weighting 
the field galaxies with the mass and star formation rate supernova frequencies of the simple A-l-B 
model produces good agreement with the observed clustering. We conclude that these supernova 
clustering differences are primarily the expected outcome of the dependence of supernova rates on 
galaxy masses and stellar populations with their clustering environment. 
Subject headings: surveys - supernovae: general - galaxies: clustering 



1. INTRODUCTION 

The relationship between clustering environment, 
supernova production rates and subsequent feedback 
to the environment has long been of interest. The 
early investigations of the relation between la super- 
nova rates and clustering were based on local samples 
largely targeted towards clusters. A study of 40 
some nearby supernovae primarily in cluster fields 
concluded that field and cluster rates were not distin- 
guishable (jZwickv 19421 : iBarbon 19681 ). Later, when 
the host galaxy luminosity dependence of supernova 
production frequencies was recognized, t he cluster rat e 
was found to be lower than the field (iBarbon 19781 ). 
Most supernovae in those samples are likely Type la 
although the data and supernova typing methods of 
the time were not able to make a clean separation 
between core-collapse supernovae and white dwarf 
explosions. The direct dependence of the supernova 
rate on luminosity, or total stellar mass and the star 
formation r ate was remarked on in a number of different 
contexts (Ivan den Bergh 19591: Ivan den Bergh 196C : 



Tammann 1970t 



ICrane Tammann fc Woltjer 1977 : 
Caldwell fc Oemler 1981t ). One recent search targeted 
more than one hundred Abell clusters found six new su- 
pernova la and concluded that clusters had a supernova 
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rate similar to the field (| Sharon et al. 20071) . A recent 
analysis of a sample of 136 local supernova reached 
the conclusion that cluster early type galaxies have a 
supernova la rate about three times the rate in fiel d 
early types at 98% significance (iMannucci et al. 2007f) . 
The application of supernovae as cosmological distance 
indicators is subject to the concern that the intrinsic 
luminosities are evolving. These clustering studies raise 
the possibility that supernovae, and especially those 
in elliptical galaxies, may be subject to environmental 
influences. This paper first undertakes a new evaluation 
of the influence of clustering and then tries to determine 
whether any differences can be reasonably explained as 
the result of host galaxy population differences with 
clustering environment, i.e. the well known tendency 
for relatively more early type galaxies to be present in 
clusters. 

The Supernova Legacy Survey (SNLS) sample probes 
out to redshift one, a range which contains substan- 
tial evolution of the host galaxy and supernova popu- 
lation which provides guidance to processes at higher 
redshifts when galaxies were undergoing their primary 
star formation. Supernova la are of particular interest 
because they dominate the production of iron peak ele- 
ments and inject some 10^^ erg of kinetic energy start- 
ing approximately 10^ years after their progenitor stars 
were formed. A fleld supernova survey like SNLS of- 
fers a fair sample over all clustering environments but 
has the statistical difficulty that only a few percent of 
all galaxies are in clusters, leading to relatively small 
numbers even in relatively large sample s like the SNLS 
(|Astier et al. 20061 : iGraham et al. 20071) . However, the 
SNLS has the density and depth of both supernova 
hosts and field galaxies that we can undertake a cross- 
correlation analysis that takes advantage of the approxi- 
mately 280 identified host galaxies in our current 3-year 
sample. This type of analysis does not separate the en- 
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vironment into field, group or cluster, but usefully ex- 
presses the statistical average over all environments and 
complements analyses which specifically look at clusters. 

A simple, testable, theory relating supernova hosts to 
field galaxies is available. Quite generally the rate of 
production of la supernovae is the convolution of the 
star formation rate with D{t — tf), the fractional dis- 
tribution of delay time between the formation of the 
stars at tf and that fraction that explode as la at t, 

or, R{t) = /q D{t - tf)M(tf)dtf. The delay time distri- 
bution can, in principle, be calculated given the evolu- 
tion and (largely u nknown) frequ encies of close binaries 
into supernova la (iGreggio 2005f) . Scannapieco & Bild- 
sten(2005) introduced the simple empirical approxima- 
tion that breaks the delay time distribution into a two 
piece discontinuous model having a prompt component 
proportional to a star formation indicator and a delayed 
component proportional to the built up stellar mass of 
the galaxy. That is, D{t-tf) is A for [0, t-At], and B/At 
for [t — At,t]. Doing the integrals gives the resulting su- 
pernova la rate, as R{Ia) = AM -\- B{M)At, where the 
stellar mass, M, and recent star formation rate, {M)At 
(abbreviated to M) can both be inferred from our data. 
The model can be fit to the host and field galaxies at a 
single redshi ft, but then predicts a redshift dependence 
of the rates (ISulliyan et al. 20061 : iMannucci et al. 20071 : 
iHowell et al. 2007D and the clustering dependence. In 
this paper we will test whether this weighting of an ap- 
propriate field sample can reproduce the clustering of 
supernova la hosts. 

2. DATA AND METHODS 

The identifications of the supernova host galaxies and 
the mass and star formation rates of all galaxies are 
derived with the methods of Sullivan et al. (2006). In 
brief, we have deep [um, <7m, ^m, ^m] images of the 
four SNLS fields in the Megaprime filter system (similar 
to the SDSS filters). From the photometry we derived 
indicative star formation rates, stellar masses and pho- 
tometric redshift s for all of the galaxies using the PE- 
GASE approach (iLe Borgne fc Rocca-Volmerange 20021 : 
iLe Borgne et al. 2004D . The SNLS sample used here in- 
cludes detections up to early 2007 with 281 spectroscop- 
ically confirmed supernova la host galaxies. To produce 
a sample suitable for clustering studies, we first ensure 
that field and hosts cover the same observational space in 
brightness, redshift, and mass and on the sky. The first 
cuts are simply to remove supernovae that have no clear 
host or are in regions of the array that are masked out, 
usually because of the presence of a bright star. Second, 
we restrict the mass of the host galaxy to be in the range 
of 10^ M© to lO^i-^^ Mq which ehminates about twenty 
low mass hosts, but none at high mass. The main bene- 
fit is to remove about half of the field sample and hence 
reduce the noise in the measurement. We restrict the 
galaxy brightness range to 20 < < 25 mag. Figure 1, 
and the redshifts to be in the range of 0.2 to 0.9, Figure 2. 
The sample starts to become incomplete around redshift 
0.6, but this is not a problem for a correlation analysis. 
The cuts leave a total of 163 supernova host galaxies and 
147,246 field galaxies in the four fields. The 46 hosts in 
the range 18 < im < 20 mag are at low redshift and com- 
plicate the problem of producing a comparable reference 
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Fig. 1. — The im magnitude distribution of host galaxies 
(soHd/red), imitation hosts (dashed/green, offset left), and field 
galaxies (dotted/blue, offset right). Each histogram area in the 0.5 
mag bins is normalized to unity. 
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Fig. 2. — The redshift distribution of host galaxies (solid/red), 
imitation hosts (dashed/green, offset left), and field galaxies (dot- 
ted/blue, offset right). Each histogram areas in the Az = 0.05 bins 
is normalized to unity. 

population from the field galaxies, which are mainly at 
higher redshifts. 

We will measure the clustering of galaxy hosts relative 
to the field with the angular two-point cross-correlation 
function. The two-point function gives the sky density 
of field galaxies around a supernova la host galaxy at an 
angular separation 9 as nhf{0) = n/[l -I- Whf{d)\- For 
galaxy magnitudes in the range of an ~ 25 mag clus- 
tering of field galaxies leads us to expect that the Wh / {0) 
will be generally less than one beyond a few arc-seconds, 
that is, over the entire range we probe. To estimate 
these weak correlations we use the cross-correlation form 
of the Landy-Szalay (1993) estimator modified to have 
the required symmetry in the two populations (Blake 
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The random sample is generated by uniformly plac- 
ing 10^ points on the sky in each field and then remov- 
ing the points that fall in the masked out regions. This 
leads to about 90,000 points per field, about three times 
the number of field galaxies, which is sufficient to re- 
duce shot noise errors well below field to field variance. 
We take advantage of the spectroscopic redshifts of the 
supernovae hosts and the photometric redshifts of the 
field to boost the signal to noise of the clustering mea- 
surement by requiring that the field galaxy be within 
±0.1 in redshift of the host spectroscopic redshift. The 
selected redshift window approxim ately matches the er - 
rors in the photometric redshifts (jSuUivan et al. 20061 ). 
Other values were tried, but the selected 0.1 is close to 
optimal, giving about a factor of two improvement over 
no redshift selection at all. 

To compare the clustering of supernovae hosts to field 
galaxies, we construct an imitation host sample dis- 
tributed in redshift and brightness identically to the real 
hosts. We randomly draw from the field galaxies a sam- 
ple of 10000 galaxies in each field with the same im dis- 
tribution as the hosts. Figure [T] The figure shows that 
the supernova hosts have a very broad flat brightness 
distribution compared to the rising counts of the field 
brightness distribution. Once the magnitude selection is 
done, the redshift distribution of the imitation hosts is 
sufficiently similar to the host galaxy redshift distribu- 
tion as shown in Figure [2] that we do not futher select 
to improve the imitation hosts. An important aspect of 
this sample is that the precision of the clustering mea- 
surement derived from it is completely dominated by the 
structures in the 4 fields. To quantify this, we reduced 
the sample to 3000 imitation hosts per field, finding that 
the errors have increased by 13%, far short of the 
expected from sample size dominated Poisson statistics. 
We conclude that for lO'' imitation hosts the Poisson er- 
rors are about a 10% component of the errors, with the 
field to field component being in common to both sam- 
ples. 

We will compare the host-field cross-correlation with 
the imitation host-field cross-correlation. The imitation 
hosts are compared with unit weights and with weights 
assigned from the A-l-B formulation which reproduces 
the expected supernov a frequ encies. The A and B val- 
ues from (Sullivan et al . 20061 ) are adopted. In principle 
these are not quite appropriate, since the hosts here are 
not weighted for detection efficiency. We did determine 
A and B for the sample used here, finding no practical 
difference in the outcome. 

3. RESULTS 

Figure [3] shows the angular cross-correlation of the 
field and hosts in the redshift range 0.2 to 0.9, with 
all galaxies required to have 20 < < 25 mag and 
8 < logAf* < 11.55. On the average, the host-field cross- 
correlation is slightly stronger than the comparably se- 
lected field galaxies, with the difference being largest at 
separations closer than about 10". To quantify the cross- 
correlation we fit the power law model w{6) = Aiu9~^ , 
equivalent to {9o/0)~^, where is the correlation an- 
gle. The fitting is done as a linear fit in logw — \og9. 
The error at any point is estimated with the Jackknife 
procedure, in which the difference of the results recalcu- 
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Fig. 3. — The log-log plot of the angular cross-correlation be- 
tween the field galaxies and the supernova hosts (dashed/red, open 
circle), the imitation hosts population (dashed-dotted/green, open 
square) and the A+B weighted imitation hosts (dotted/blue, dia- 
monds). 
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lated dropping one field at a time and all four fields are 
summed in quadrature and divided by n — 1 to estimate 
the variance. The jackknife approach gives realistic, un- 
biased errors. We exclude measurements inside 3" on the 
basis that the galaxies begin to physically overlap and it 
removes the self-correlation component. We fit over the 
3— 100" range where the correlations are always positive. 
On the average the correlation of the real hosts is 60% 
stronger over the entire radial range, a full factor of three 
more in the 3— 10" range about about 12% stronger over 
the 10 - 100" range. 

Table 1 gives the best fit [log A^, , S] values for the 
power-law fit to 3 — 100" range of the angular two-point 
correlation functions in the first two columns. The fit- 
ted parameters are shown in Figure 4, along with their 
error ellipses, which show the usual strong correlation 
between the amplitude and slope. The third column of 
Table 1 gives the reduced of the cross-correlations 
when modelled with the host field galaxy fit in the first 
row. The true host fit is rejected as inconsistent with the 
imitation hosts at 99.3% confidence. A better measure 
of the significance of the clustering difference between 
the hosts and the field is the difference of the two 
correlations functions. As discussed above we calculate 
using the jackknife computed variance of the host- 
field cross-correlations which includes both the Poisson 
errors and field-to-field variance which is in common, as 
discussed above. These values are presented in the last 
column of Table 1. We see that the hosts and imitation 
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Fig. 4. — The fit of the measured correlations to = Sg^ with 8 
measured in arc-seconds. Power law fits over the range 3 — 100" are 
shown for the supernova hosts (red, open circle), imitation hosts 
(green, open square), and A+B weighted imitation hosts (blue, 
diamond), A=0 (cyan, filled circle) and B=0 magenta, hexagon). 
The error ellipse is computed from the co-variance matrix of tfie 
fit. 

hosts have a per 8 degrees of freedom of 4.3, which has 
a probabihty of occurrence less tha n 3 x 10~^. The ta- 
ble also presents the results for (Sullivan et al. 2006l )the 
A-l-B weighting of the host galaxies in the correlation 
function. Weighting with star formation or mass alone, 
A — 0, partially explains the effect, but the reduced of 
the difference between the correlation functions indicates 
that the agreement is not statistically adequate. Mass 
weighting alone does provide an acceptable fit to the clus- 
tering pattern. We conclude that the A-l-B weights work 
to account for the rate of supernovae in strongly clus- 
tered, relatively high mass, galaxies. 

To estimate the physical correlation length we use 
the Limber equation. The photometric redshift dis- 
tribution is adequately fit with n{z) proportional to 
0exp(-i[(z-z)/(72]2) with z = 0.6 and cr^ =0.18. This 
distribution is convolved with a Gaussian with a redshift 
dispersion of 0.1. We find that the host-field correlation 
angle of approximately 1.5 ± 0.5" with S — —0.63 ± 0.13 
corresponds to an (assumed constant) co-moving correla- 
tion length of 6.5zt2h~^ Mpc in a flat ^Im = 0.24 cosmo- 
logical model. It must be borne in mind that this cross- 
correlation applies to a specific field sample, spanning the 
10^ — 10^^-^^ M0 mass range. The correlation length is 



marginally stronger than the canonical 5h~^ Mpc of typ- 
ical luminosity galaxies, and bolsters the view that the 
relatively more massive galaxies dominate the supernova 
clustering. 

4. DISCUSSION 

The main outcome of this paper is the first measure- 
ment of the two-point correlation between supernova 
hosts and the surrounding field galaxies. Over this red- 
shift range we find that the hosts are more correlated 
than a similar zm magnitude and redshift selected field 
sample with the probability that this is a chance event 
of 3 X lO"-"^. Weighting the field with AM^ + BM re- 
produces the host clustering remarkably well, primarily 
as a result of accounting for the higher mass-to-light of 
older stellar populations. Mass weighting alone can ex- 
plain the clustering properties for this sample as shown 
in Table 1, although it is important to note that the 
supernova rates at these redshifts are completely domi- 
nated by the star formation rates, not the mass, of the 
galaxies. The star formation rate is correlated with the 
galaxy mass which partially explains why either mass or 
star formation weighting alone works, but the very strong 
correlations of the most massive galaxies, which have ef- 
fectively no star formation but a significant supernova la 
production, is also part of the effect. We conclude that 
supernova la host galaxy clustering is primarily the out- 
come of the dependence of galaxy stellar populations on 
clustering environment. 
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